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Abstract. The vacuum (Casimir) energy in quantum field theory is a problem relevant both to new 
nanotechnology devices and to dark energy in cosmology. The crucial question is the dependence 
of the energy on the system geometry under study. Despite much progress since the first prediction 
of the Casimir effect in 1948 and its subsequent experimental verification in simple geometries, 
even the sign of the force in nontrivial situations is still a matter of controversy. Mathematically, 
vacuum energy fits squarely into the spectral theory of second-order self-adjoint elliptic linear 
differential operators. Specifically, one promising approach is based on the small-f asymptotics 
of the cylinder kernel e^'^, where H is the self-adjoint operator under study. In contrast with the 
well-studied heat kernel e^'^ , the cylinder kernel depends in a non-local way on the geometry of 
the problem. We discuss some results by the Louisiana-Oklahoma-Texas collaboration on vacuum 
energy in model systems, including quantum graphs and two-dimensional cavities. The results may 
shed light on general questions, including the relationship between vacuum energy and periodic or 
closed classical orbits, and the contribution to vacuum energy of boundaries, edges, and corners. 

Keywords: Casimir force, vacuum energy, quantum graphs, quantum billiards 
PACS: 03.70.H-k, ll.lO.Gh, ll.80.La, 03 .65.Sq 

1. INTRODUCTION 

Since Casimir's famous calculation in 1948 showing an attractive force between parallel 
conducting plates due to vacuum fluctuations of the electromagnetic field [1], forces as- 
sociated with vacuum energy in quantum field theory have been studied in a wide variety 
of contexts [2]. These range from the bag model of the nucleon [3], to cosmology [4], 
to stabilization of brane world models [5], and to practical applications in micro- and 
nano-electromechanical systems [6]. 

Of course, the calculational details in specific applications will depend on the system 
dimension, the nature of the relevant quantum fluctuating field (e.g., a vector electro- 
magnetic field), and the detailed boundary conditions (e.g., ones that properly take into 
account the finite plasma frequency in the electromagnetic case). In the examples consid- 
ered here, we bypass these application- specific details and instead consider a toy model 
of a scalar field, usually with Dirichlet or Neumann boundary conditions. As we will 
see, these simple examples will help to elucidate important general questions concern- 
ing Casimir forces that are independent of the specific context. These questions relate 
to proper regularization and renormalization of the formally infinite vacuum energy, the 
relation of Casimir forces to periodic and closed classical paths, and the role of bound- 
aries, edges, and comers. 

Formally, the vacuum energy of a scalar field is given by §Lnft)„, where (0„ are 



the eigenfrequencies, given by solutions of — V^9„ = ^(p„ with the relevant boundary 
conditions. In the following, we work in units where h = c = I. Of multiple methods of 
regularizing the infinite vacuum energy (including e.g., dimensional regularization), we 
focus here on the time-splitting regulator, associated with the cylinder kernel Tt{x,y) = 

{x\ e^^~'^^' \y) . Defining 

the physical vacuum energy is given by taking the limit of as ? — >■ if this limit exists 
(i.e. if the divergent terms can be shown to cancel). 

2. VACUUM ENERGY IN QUANTUM GRAPHS 

We begin by applying the above approach to quantum graphs, a class of one-dimensional 
models that have been widely used as approximations for the free-electron theory of 
conjugated molecules in chemistry, for quantum wire circuits in nanotechnology, and 
for photonic crystals in optics. More generally, quantum graphs provide a useful testing 
ground for investigating general questions about quantum chaos, scattering, and spectral 
theory. A good discussion may be found in a recent survey by Kuchment [7]. 

Mathematically, a quantum graph consists of one-dimensional bonds meeting at ver- 
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tices, with the scalar field satisfying — V^^„ = on each bond and prescribed bound- 
ary conditions at each vertex. For detailed presentations of the mathematical model, see 
Refs. [8, 9]. Vacuum energy in quantum graphs has been studied recently by Berkolaiko, 
Harrison, and Wilson [10]; here we show some results obtained by Fulling, Kaplan, and 
Wilson [11]. 



2.1. Pistons in One Dimension 

The left panel of Figure 1 shows a simple line graph consisting of three bonds and 
four vertices. The two middle vertices are to be thought of as movable pistons, i.e. 
their position may change keeping the total length Li+a + L2 fixed. The objective is 
to calculate the vacuum energy as a function of these positions, and thus to obtain a 
Casimir force acting on the pistons [12]. Focusing first on the bond of length a separating 
the pistons, we note that the eigenfrequencies are (On = mt/a, where n= 1,2, •• • for 
Dirichlet boundary conditions at the pistons. We then have 
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FIGURE 1. (Left) A line graph consisting of three bonds and four vertices. The two middle vertices are 
movable pistons. (Right) A star graph with B — S pistons. 



and the regularized vacuum energy is 

^• = 7^2-i-,^0(,). (3) 
Adding the vacuum energy from the other two segments, we find 

2nt^ 24a 24Li 24Li ^ ' ^ ^ 

Thus the divergent term corresponds to a geometry-independent constant energy density, 
and is unobservable since it will not contribute to a force on the piston. After safely 
discarding this constant energy shift, we may let ? — )• and Li,2 — > °° and obtain the 
well-known finite, attractive force 
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The same result is obtained if Neumann boundary conditions obtain at each vertex 
{(On = nil I a with n = 0, 1,2,---)- However, if one piston is Dirichlet and the other 
Neumann, the analogous calculation yields a repulsive force 

F„ = + j|^. (6) 

Though simple, the calculations yield little or no insight as to why the force may be 
attractive in some situations and repulsive in others. To obtain such insight, we turn to 
an alternative perspective. We first note that 

TrTf = JdxTt(x,x). (7) 



Now the free cylinder kernel in one dimension is 

T,%x,y) = -- \. ^. (8) 

Then Tt (x, x) in a problem with boundaries is obtainable by the method of images as a 
sum over periodic and closed orbits: 

/ A 

Tt (x, x) = Re y — + closed orbits , (9) 

where p labels periodic orbits passing through the orbit length, and Ap is the 

product of scattering factors at vertices (e.g., —1 at a Dirichlet vertex and +1 at a 
Neumann vertex.) The expression arising from closed orbits (i.e. orbits starting and 
ending at x but with opposite momenta) is similar, and is omitted here because after 
integration over x, closed orbits in graphs may be shown to give zero contribution to the 
total energy. This is not the case in two- or three-dimensional billiards, as discussed in 
Section 3. 

In our case, all periodic orbits in the region between the pistons are repetitions of the 
primitive periodic orbit of length 2a. Separating out the zero-length orbit (r = 0) and 
taking the trace, we find 

(10) 



Trr, = |jxr,(x,x) = ^^+Re£^^ 



ra)2 + ?2 

where r labels the repetition number and A is the product of scattering factors for the 
primitive two-bounce orbit. 

The asymptotic t ^0 behavior may now be analyzed term by term. We immediately 
see that all orbits of nonzero length make contributions ~ ? to the cylinder kernel, and 
thus finite independent contributions to the energy. The t ^ divergence is seen to 
be associated exclusively with the zero-length periodic orbits, which exist at every point 
X and yield the divergent but constant and geometry-independent energy density. 

The periodic orbit sum converges. Differentiating Tr Tj with respect to t to obtain the 
vacuum energy and then with respect to a to obtain the force on a piston, we have 

1 °° A'' 

We note that A = + 1 if the pistons are both Neumann or both Dirichlet, and A = — 1 
for mixed boundary conditions. Eq. (11) thus reproduces the results of Eqs. (5) and (6); 
moreover we see that the sign of the force arises from reflection phases, and is already 
correctly obtained if we consider only the phase associated with the shortest periodic 
orbit (r = 1). 



2.2. Pistons in General Star Graphs 



We turn our attention to star graphs, an example of which is pictured in the right panel 
of Fig. 1. A total of B line segments of large length L meet at the central vertex, where 



the wave function satisfies Kirchhoff boundary conditions: (i) continuity (pj{0) = ^/t(0) 
for all j.,k — l - B and (ii) current conservation JL^^i (Pj{0) — 0, where ^j(O) is the 
outward derivative along bond j. Along each segment j, a piston is located at distance 
Gj from the central vertex, and the boundary condition imposed by the piston may be 
Dirichlet (reflection with phase —1), Neumann (reflection with phase +1), or reflection 
with an arbitrary phase e'^j (to break time reversal symmetry). We will be interested in 
computing the dependence of the energy on the piston positions, i.e. in the Casimir force 
on the pistons. 

We focus initially on the "interior" of our system, i.e. on the graph consisting of B 
bonds of length ai ■ --aB, and excluding the space beyond the pistons. For B > 2 and 
generic aj, no analytic expression exists for the spectrum or for the vacuum energy, and 
a numerical approach must be employed. For a general quantum graph, the eigenfre- 
quencies are given by solutions of a characteristic equation det /z((o) = [8]; in the case 
of a star graph with irrationally related bond lengths aj the relevant equation becomes 

B 

J^tan((jOa^ + 0y) = 0, (12) 
;=i 

where dj = or ;r for a Neumann or Dirichlet piston on bond j, respectively [11]. If 
we numerically obtain in this way all eigenfrequencies (On up to a cutoff (Bbax, we may 
write 

£™°i=l £ (One-"^' = Et + 0{e-'^'^'). (13) 

Since the error associated with omitting eigenfrequencies (0„ > OVax is 0(e^®™'''^), we 
must consider (Omaxt ^ 1. 

Now we turn to the "outside," i.e. the segments aj <Xj <L located beyond the pistons. 
From Section 2.1 we know that the outside energy consists of finite terms that decay as 
1/L and may therefore be neglected for large L, plus a divergent l/t^ term associated 
with a geometry-independent constant energy density. The divergent terms, as before, 
will combine with the divergent part of the interior vacuum energy to yield a constant 
energy shift BL/2nt^ proportional to the total length BL and independent of the piston 
positions. To obtain the physical forces on the pistons for large L we thus need only to 
subtract from the interior energy the divergence proportional to the total interior length 
X)y flj. The physically observable energy is then given by 

£fnite^^num_^Weyl^ ^^^^ 

where the divergence coming from integrating the Weyl density in one dimension 
p{(o) — YijCti/'^ between and (Omax is 



Expressing the finite part of Et as a power series. 



Ef'''^ = EQ + ait + a2t^ + - 



(16) 




FIGURE 2. (Left) The force on a piston in a star graph with B bonds of length 1 and either Neumann 
or Dirichlet boundary conditions at each piston is computed using only repetitions of the shortest periodic 
orbit and compared with the exact answer. Positive values indicate repulsive forces. (Right) The difference 
between the exact Casimir energy Eq and a periodic orbit sum including all orbits up to length Lmax is 
shown for a star graph with four bonds of length 1.1, 1.6176, 1.2985, and 1.1159, and a Neumann piston 
at the end of each bond [11]. 



and numeiically evaluating Ef for several values of the cutoff t with CO^^y^ <^ t <^ 
min(aj), we easily obtain the vacuum energy Eq for any given star graph to any desired 
order of accuracy. 

To attain a more physical understanding, we wish to compare the numerical results 
with periodic orbit sums. For example, taking all pistons to have Neumann boundary 
conditions and summing only over repetitions of the shortest periodic orbits (i.e., over 
orbits bouncing back and forth in a single bond), we obtain 

pshon..,_ " {, 241n2 \ » 1 

which for large B compares well to the analytic result — g)ffor5 equal-length 
bonds. Similar results are obtained in the Dirichlet case. As illustrated in the left panel of 
Figure 2, both repulsive behavior in the Neumann case and the attractive behavior in the 
Dirichlet case are well explained by considering only the shortest periodic orbit. For a 
better quantitative approximation we may add contributions from orbits of longer length. 
Figure 2 (right panel) shows the convergence to the numerically exact Caimir energy Eq 
for a typical star graph with Neumann pistols. Here the error scales as L^l^^, where Lmax 
is the length of the longest orbit included. For mixed Dirichlet or Neummann boundaries, 

—3 /2 

or for arbitrary phases at the pistons, the convergence is shown to be Lmax [1 !]• 



3. VACUUM ENERGY IN RECTANGLES, PISTONS, AND 

PISTOLS 

We now extend the approach of Section 2 to two-dimensional billiards (the extension to 

the three-dimensional case and the electromagnetic field is also straightforward [13].) 
An important motivation for this work [14] is to investigate the physical reality of the 



outward force on the walls of a square or cubic box, as obtained by Lukozs using naive 
renormalization (i.e., by simply discarding infinite terms) and ignoring the outside of the 
box [15]. 

We begin with a rectangle of sides a and b. As for a line segment (Section 2.1), we 
can use the method of images to evaluate Tr Tt, and thus the regularized vacuum energy 
Ef. Each term in the image sum may be associated with a classical path leading from x 
to X in the rectangular, and these terms may be classified by the number of bounces the 
path makes off the vertical and horizontal walls. Periodic paths make an even number of 
bounces 2 j off the vertical sides and an even number of bounces off the horizontal sides. 
The resulting contribution to the vacuum energy is 



^ [?2+(27fl)2+(2)t&)2]5/2' ^ ^ 

where r]j]^ is the number of Dirichlet bounces for a given orbit, and we have separated out 
the purely vertical and purely horizontal periodic orbits (7 = and k = 0, respectively), 
as well as the zero-length orbit 7 = = 0. As in the one-dimensional case, the zero-length 
orbit contributes a divergent but constant and geometry-independent energy density, i.e., 
a divergence proportional to the total area ab. Assuming all Neumann or all Dirichlet 
sides, so that all r\j]^ = 0, we have 

ab ( a b\ ab ^ ^ -2 , , 2, 2N-3/2 . 2n ..q. 

^.,Periodic=2^-^(^^ + ^J-g^LL(a7 +^^) +0{t). (19) 

In contrast with the one-dimensional case, here the non-periodic closed orbits (ones that 
make an odd number of bounces off the horizontal sides, or an odd number of bounces 
off the vertical sides, or both), make a nonzero contribution to the total energy, including 
a divergent contribution proportional to the system perimeter. Combining periodic and 
non-periodic terms we obtain 

Area Perimeter ^(3) / a b 



27tt^ Sm^ \6% \b^ ^2 
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j,k=i \ / 

where =f refers to Neumann or Dirichlet boundaries, respectively. Naively discarding 
the divergent terms and differentiating with respect to a we find an attractive force for 
a<^b (as expected in the limit of two infinite parallel plates), but a repulsive force for 
the square a = b. 

The above analysis has two (related) problems: discarding divergent terms and ignor- 
ing the outside of the box. Both may be cured by considering a piston configuration [12], 
as discussed in Section 2 and illustrated in Figure 3 (Left). Adding contributions from 
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FIGURE 3. (Left) A piston in a rectangular box. (Right) A pistol configuration, consisting of a "barrel" 
and a "bullet." 



the a X rectangle and the {L—a)xb rectangle, we see that the divergent terms combine 
to yield contributions proportional to the total system area, and total boundary length, 
and thus independent of the position of the piston. Other than these divergent terms, 
the only contribution from the (L — a) x b rectangle that survives the L — )■ oo limit is 

(^F^) ' combining this term with the finite part of Eq. (20) and differentiating with 
respect to a we obtain a finite Casimir force, 

^piston = ^ E {2Kjk^^ . (21) 

This force is always attractive (decaying exponentially for a b and reducing to the 
parallel plate limit for a <S b). 

Finally, to approach the original motivating situation of a box with a loose lid [15] 
and to address the question of what happens when the external shaft is not present, we 
consider the "pistol" configuration illustrated in the Right panel of Figure 3. Here all 
system dimensions other than possibly the gap c are assumed to be large compared to 
the ultraviolet cutoff t. We then use scaled variables c ~rt,a = st,b = ut, d= {£ — s)t, 
where s,u,£:^ 1, and for all Dirichlet boundaries obtain 



_ us y l-2k^U^ us y I-27V 

P'^'°i ~ (1 +4k^u^y/2 '^mp^ (1 +47V)5/2 

^ lit p (l+47V + 4;t2„2)5/2 

s y -1+4;V lr{l-s) y \-lk^r^ 

In the case of a narrow chamber, a <^ b^l^c^/^ , we recover an attracive force ~ 1/a^, as 
required in the parallel plate limit. In the opposite case of a long chamber, a ^ b^/^c^/^ , 
we find that the gaps of width c dominate and we obtain an a-independent force that is 
attractive for c > at and repulsive for x < at, where a ^ 0.5888. This last situation, 
however, is least convincing physically, as we need to be in a regime where the gap 
dimension is smaller than the cutoff. See Ref. [14] for a detailed discussion. 



4. GENERAL BILLIARDS 



We note that the numerical approach to calculating the vacuum self-energy, applied to 
general quantum star graphs in Section 2.2, may be equally well applied to two- or three- 
dimensional systems, provided the spectrum may accurately be computed numerically 
up to some maximum frequency CtVax- Of course the appropriate Weyl term containing 
a\\ t ^ divergences must be subtracted from the numerical sum (13) before the 
numerical limit t <^ I may be considered. For example in the case of the interior 
of a three-dimensional cylinder with polygonal cross sections and Dirichlet boundary 
conditions [16], Eq. (15) becomes 



' 2j-co2nJo \ 2n 4n J 



1 f°° dk /"^ax fTTTTT?. -tVl?+a^ ( /Area Perimeter^ 

(23) 

where the area and the perimeter refer to the polygonal cross section, C = Y,i ~ 

with a,- the interior comer angles of the polygon, and (Umax is the maximal eigenvalue 
obtained numerically for the polygon. An additional divergent term In? must be consid- 
ered in the presence of boundary curvature. These approaches are now being applied to 
study the self-energy of stadium-shaped and elliptical cavities, as well as to investigate 
the self-energy associated with the outside of a billiard of arbitrary shape. 



5. CONCLUSIONS 

We have seen that careful regularization and renormalization (including both inside and 
outside contributions) are needed to obtain physically meaningful vacuum energies and 
Casimir forces. Classical orbit approaches, including both periodic and non-peridoic 
orbits, produce exact results in simple cases and may allow for good approximations 
where exact solutions are nonexistent, including general quantum graphs and polygonal 
billiards. Furthermore, such semiclassical approximations may be compared with results 
obtained numerically by directly summing eigenfrequencies and subtracting known 
divergences associated with the Weyl part of the spectrum. Intelligent combination of 
analytical and numerical tools can be a promising tool for furthering our understanding 
of Casimir forces in general geometries. 
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